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Abstract 

A finite number of gravitating point particles in 2+1 dimensions may close the universe 
they are in. A formalism previously introduced by the author using tesselated Cauchy 
surfaces is applied to define a quantized version of this model. Special emphasis is put on 
unitarity and uniqueness of the evolution operator and on the physical interpretation of 
the model. As far as we know this is the first model whose quantum version automatically 
discretizes time. But also spacelike distances are discretized in a very special way. 



1. Introduction 

Consider the Lagrangian 

£ = ^V=gR + V^gi-b^W*** - X^ 2 ) ( L1 ) 

in 2+1 dimensions. The perturbation expansion in Newton's constant G is non- 
renormalizable. Yet one may suspect that a quantum version of this model exists, because 
in a very special classical limit the system is integrable. We have in mind the limit where 
the <p particles with mass m become classical gravitating point particles [1] . There are no 
gravitons in 2+1 dimensions [2,3]. The particles move in rectilinear orbits in a locally flat 
space, and this motion is non-trivial only because the continuation of these orbits depends 
on an element of a braid group. 

But the small distance behavior of a quantum field theory described by (1.1) must 
be very peculiar. Newton's constant G defines a Planck length, and at this length scale 
any perturbative approach will break down. Typical quantum gravitational effects must 
be essential there. 

Now the pure gravity system, without particles, can be quantized [4]. But adding 
spinless point particles is essentially equivalent to adding a scalar field, and this may 
provide us with an infinite dimensional Fock space. If second quantization should be 
inevitable we should prepare for creation and annihilation of particles, and this was not 
considered in [4]. 

In previous papers [5] we expressed doubt whether a rigorous quantum version exists 
at all, because of the requirement of a fundamental quantized gravity theory. But further 
examination of the classical system provided us with so much information concerning its 
fundamental degrees of freedom and its causal structure that a renewed attack is made 
possible. We here report on a beautiful short distance structure, showing features of 
finiteness that were previously speculated upon for the 4 dimensional theory. 

Our approach for the classical system will be that of Ref. [3] . Let us briefly summarise 
the method. We start with a set of Cauchy surfaces C . A Cauchy surface is a spacelike 
cross section of space-time. Here it is 2 dimensional. No pair of points in a Cauchy surface 
is allowed to be timelike separated, and all other points in space-time must be either in the 
future or in the past of at least one point in each Cauchy surface. To each of our Cauchy 
surfaces C we attach a time parameter tc ■ 

Next, we design a tesselation of each Cauchy surface, so that its evolution to the 
future (and to the past) can be calculated. Each of these "tiles" is a polygon. Particles 
may only sit at the corners of the polygons, so inside each polygon space-time is flat. The 
constant time surface t = tc defines a preferred Lorentz frame (but not yet its origin) 
for space-time there. The polygons are bounded by edges. At an edge, the two Lorentz 
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frames of the adjacent polygons are related by a Lorentz transformation. Since at an edge 
the two time coordinates coincide the Lorentz boost from one polygon into the next must 
be directed orthogonally to this edge. The boost parameter for an edge Lj will be called 
2rji. For reasons that will become clear later we will now choose the signs such that rji > 
if both polygons contract; rji < if they expand. The velocity of the edge itself in both 
frames is described by half this boost, rji . A particle at a corner of a polygon will connect 
two edges that will be glued together in such a way that the particle produces a conical 
singularity. In general the particle moves, so that the Lorentz frames at both sides of such 
an edge will also be related by a Lorentz transformation. If 2/3 is the deficit angle at a 
particle, m is its mass, and £ the boost parameter for its velocity in the local frame, then 
we have [3] 

tan/3 = cosh£ tan |m , 
tanhry = sin/3 tanh£ , 
cos \m = cos (3 cosh rj , 
sinh?] = sin |m sinh£ . 

(for future calculations it turned out to be convenient to absorb here a factor 2n in the 
definition of the masses m of the particles, as compared to our earlier expressions in ref 
[2-5]). 

The topological structure of a tesselation will be denoted by a diagram indicating 
the edges of the polygons without bothering about actual lengths or angles. Depending 
on the global topology of 2-space the diagram should be seen as living on a topologically 
non-trivial sheet, which we unfold by removing a few points. The diagram (after adding 
the point(s) at infinity) indicates how the polygons fit together. 

The evolution is now indicated diagrammatically. During short intervals of time we 
may simply allow time to evolve equally fast on all polygons, so that the edges move with 
their well-defined velocities. But it will be unavoidable that as time continues something 
will happen. It could be that the length of an edge shrinks to zero. It could also happen, 
since many polygons are not convex, that one of the vertices of a polygon hits one of the 
other edges, at which point also it becomes illegal to continue the description in terms of 
these particular polygons. A transition in terms of another set of polygons takes place. It 
is the succession of many such transition that we studied. The complete set of all possible 
transitions in a diagram is listed in Fig. 1. 

In most cases a new edge is created, which implies that two polygons that were truly 
separated before, now will acquire an edge in common, whereas other edges may disappear. 
Since the relative Lorentz transformation between one polygon an an adjacent one was 
determined by the succession of Lorentz transformations at other edges, and since this will 
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Fig. 1. The nine different possible transitions diagrammatically 

not change, one will always be able to compute both the orientation of the new edge L\ 
relative to the others and the new Lorentz boost parameter r\\. 

In practice we compute these new numbers using "triangle relations". Consider a 
vertex between three polygons, /, II and III, and let 0^1,2,3 be the angles between two 
edges in each polygon, and 771,2,3 the three Lorentz boosts, labled as shown in Fig. 2. We 
define 



sincti = Si , coscti = q , sinh2?yj = Oi , cosh 2^ = 7^ ; (1-3) 
then we have the relations 

si : s 2 : S3 = : a 2 : 03 , (1.4) 

72«3 + S1C2 + C1S273 = , (1.5) 

ci = c 2 c 3 - 71S2S3 , (1.6) 

7i = 7273 + 0203C1 , (1.7) 

cota 2 = — cotcti cosh2?73 — coth2?72 sinh2r7 3 / sincti , (1.8) 

and all cyclic permutations. 



The use of these relations is described in detail in Ref. [3] . Any possible ambiguity in 
the parameters of a newly opened edge is removed by requiring that the edge grows with 
a positive time derivative and that the complete set of polygons must form a true Cauchy 
surface at all times. An edge L grows or shrinks at its two end points A and B : 

L = g A + g B , (1.9) 
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Fig. 2. The angles cti and the boosts rji at a vertex between three polygons 
At a vertex A the growth gA of edge L\ is given by 

#a,i = (vi cosa 3 + u 2 )/ sina 3 = (vi cosa 2 + U3)/ sina 2 , (1-10) 

where 

= -tanh//i = -cri/(l + 74) . (1-11) 
At a particle P the contribution to the time dependence L is 

<7p = tanh?] cot/5 = tanh£ cos/3 . (1-12) 

The equations (1.11) and (1.12) show how the edges evolve. 

The degrees of freedom of the system are essentially the collection of lengths Li of 
the edges % and the Lorentz boost parameters rji at these edges. The orientations of these 
edges, and with them the orientation of the Lorentz boosts there, are then fixed because 
one can compute the angles cti using first (1.7), after which any ambiguity for the sign of 
s\ can be lifted using (1.4) together with the information that at each vertex at most one 
of the Si is allowed to be negative. 

There will however be constraints. Each polygon must close exactly, which implies 
that the angles at its N corners must obey 

N 

^(tt - at) = 2n , 

i=l 

(counting the contribution of a particle P as ctp = 2n 
sum of all edges must coincide with the origin: 

N 

$>,e^=0, (1.14) 

i=i 

where u>i is the orientation of the edge Lj in the frame of the polygon, to be computed 
from the angles aj . So each polygon produces three constraints altogether. 



(1.13) 

2/3). Furthermore the vectorial 
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2. Brackets 

We have the complete set of degrees of freedom, their equations of motion (1.9) - 
(1.12), and the constraints (1.13) and (1.14) on them (which are automatically preserved 
by the equations of motion). Naturally, if we wish to find a quantum version of this 
model we have to find a Hamiltonian and Poisson brackets that generate these equations 
of motion. The following construction was discovered by first studying the weak gravity 
limit, at which space-time becomes completely flat, and the particles form a Fock space 
with known expressions for energies, momenta and Poisson brackets. In this limit the 
polygons form diagrams such that particles and clusters of particles are each connected 
with lines that are oriented in such a way that they are all parallel to the total momentum 
they carry. We read off from (1.2a) that the deficit angle 2(3 for a particle coincides precisely 
with its energy in this limit (since m is infinitesimal), and from (1.2b) we see that then 2rj 
precisely corresponds to its momentum. 

A complication, of course, is that we are dealing with a cosmology, and this implies that 
the total Hamiltonian will have a fixed value. It is natural now to take as an Hamiltonian 
the combined deficit angles. More precisely, the total energy enclosed inside any closed 
contour C is the deficit angle obtained when we parallel transport the local coordinate 
frame along this curve. 

If this is taken to be the Hamiltonian then we can deduce the canonical variable 
conjugated to the length Lj of an edge i by requiring 

Li = {H, . (2.1) 

We know that this variable must be a function of the hj , the boost parameters of all edges. 
In the weak gravity limit the variable canonically conjugated to Li simply turned out to 
be 2rji . In principle one could have expected a more complicated function of the rjj in 
the strong gravity case. Suppose now that the variable conjugated to Lj is some function 
Pi{(Vj)} • Let us then compute the Poisson bracket (2.1). 

A particle P contributes to the Hamiltonian Hp = 2/3 . Therefore it contributes to 
the time derivatives Li as follows: 

bU = {H Pl L^ = dH P /d Pl = 2(dp/di 1 )(di 1 /dp i ) . (2.2) 

Eq (1.2c) gives the relation between (3, r\ and m , from which it follows that 

d(3 sinhry cosm 

— — = ~ — — = tanhry cot (3 . (2.3) 

or) cosh rj sm(3 

We see that this is exactly the velocity gp derived earlier, eq. (1.12). Hence the contribution 
of a particle to the time derivative of its cusp agrees with the Poisson bracket only if the 
variable pi canonically associated to Li is exactly 2r\i . 
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Our scheme will only be self consistent if also the contribution of the vertices as given 
in eq. (1.10) agrees with the Poisson bracket (2.1). We expect that the contribution to the 
Hamiltonian from a vertex is 

B v = H x + H 2 + H 3 + 2ix ■ Hi = -a % ; 

7i - 7j7fc (2.4) 
cos(-Hi) = — , i,j,k = 1,2,3. 

(JjCFk 

where we used the triangle relation (1.7). 
This will give 

{H u Li} = -^dai/drji = 



{H 2: L 1 } = -\da 2 /dr]i = 



Sicr 2 a 3 
73 - 7i72 

~UU. 2 1UI\\ = 

Siaicr 2 cr 3 



ru t ^ la i a 72 -7i 73 (2.5) 



s 1 a 1 a 2 (Js 

3 



^ S2cr 3 (l+7i) 



Indeed, substituting eq. (1.7) for a 2 or ct 3 in eq. (1.10) for the velocity of the edge Li at 
the vertex V we find 

1 / c 2 ai a 3 x _ 7i + 1 ~ 72 ~ 73 (2 Q) 

,l s 2 l l + 7i l+73 j s 2 a 3 (l+ 7 i) ' 

We conclude from eqs. (2.3), (2.5) and (2.6) that indeed the variables 2r]i are the 
canonical conjugates of the Lf. 

{2r H ,Lj} = 8 ij . (2.7) 
The Hamiltonian is the sum of the deficit angles, as given by (2.4) for the vertices. 

3. Constraints 

At every polygon we need to impose the condition that all angles add up to 2n as 
given by the constraint (1.13). The angles of the polygon contribute to the Hamiltonian. 
Apparently we have 

JV 

H i = 27r (! -N) = fixed , (3.1) 

i=i 

where i lables the N corners of the polygon. 

The physical interpretation of this constraint is not difficult to see. Inside each polygon 
we had been free to choose the Lorentz frame, and in particular the time coordinate for 



6 



this local frame. If we allow this polygon to evolve all by itself it is governed by this part 
of the Hamiltonian. The constraint (3.1) tells us that this is an invariance of the state of 
the system. 

The complex constraint (1.14) must correspond to invariance with respect to Lorentz 
transformations of the frame inside the polygon. The effect a Lorentz transformation inside 
one polygon has on the surrounding Li is rather complicated, so we have not checked 
explicitly whether the change generated by this constraint indeed matches this. 

Besides these there are however more subtle constraints, in the form of inequalities: 

Li > , (3.2) 

(obviously the lengths of the edges must be greater than or equal to zero). From (1.7) one 
can also deduce that the rji must satisfy a triangle inequality: 

\Vi\ < \Vj\ + \Vk\ , i, j, k = 1,2,3 . (3.3) 
We here anticipate that another set of canonical variables will be useful: 

Xi = LiSgn(i]i) , pi = \2r)i\ , (3.4) 

which, at least classically, also obey the Poisson brackets 

{Pi,xj} = 5ij . (3.5) 

Note that the question whether the commutator analogues of (2.7) and (3.5) are also 
equivalent will be not so obvious. 

The advantage of these variables will be that there is now no further constraint on 
(the sign of) Xi , whereas pi are now limited to be positive. In addition the pi satisfy the 
ordinary triangle inequality, (3.3) without the absolute value signs. 



Fig. 3. Acceptable self-overlapping polygon 

There still is another inequality which is more difficult to write down explicitly. This 
is the requirement that all pogygons must be true polygons. The point here is that the 
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angles of a polygon need not be convex. In particular a particle with mass less than n 
will give a cusp corresponding to an angle larger than ir for the polygon. Consequently 
we must be careful that the polygon does not intersect with itself, otherwise our surface 
ceases to be a genuine Cauchy surface. In the classical case, if a polygon tends to disect 
itself it will undergo a transition of the type depicted in Fig. lc or Id. 

Occasionally however we will have a polygon that does overlap with itself; yet it may 
still be an acceptable polygon if it is still a boundary of a two-dimensional space, see Fig. 3. 
It only overlaps with itself if we insist on choosing a rectangular coordinate grid within. 

4. Quantization 

In this section the commutators will be replaced by —hi times the Poisson brackets. 
We now claim that the variables Xi and pi are to be preferred rather than the Li and 2rji 
as canonical variables. This is because at Li = a transition of the form of Fig. la, b, 
e, . . . , j has to take place, which corresponds to adding a boundary condition at Lj = 0. 
When a new edge opens up it is not a priori clear how to resolve the sign ambiguity 
in the determination of the new 77 variable, and in association with that the signs of the 
trigonometric functions sine and cosine of the new angles. If we replace the Lj by quantities 
Xi we can use the signs of Xi to give these lines an orientation inwards or outwards of the 
vertex. In particular at the transitions of Fig. le and f the newly resulting wave function 
will be continuous in x (note that in these two figures the signs of 77 on the horizontal line 
segments flip). The advantage of keeping pi positive will become apparent in Section 5. 

The transitions as pictured in Fig. 1 should be seen as boundary conditions on the 
wave function. But we can also view them as providing identities for the wave function on 
different diagrams. If the wave function is known on any particular diagram we can derive 
it on any other diagram by using these identities. For instance if in Fig. la the edge shown 
at the left has a length L > then the diagram at the right may be seen as an analytic 
continuation of it, such that the new length parameter L < 0. Figs, lc, d and lg — h show 
how polygons can be added to or removed from a diagram. 

Fig. 1 was originally intended to list the transitions for a classical theory, not directly to 
formulate the quantum system. For that it will probably be more convenient to reformulate 
the rules slightly. Technically, the transitions c and d can be obtained most easily by first 
aplitting a polygon in two, using a new edge with 77 = 0. Since then also p = this simply 
means that the wave function does not depend on the new x variable at all, and so this is 
a dummy variable at this stage. But then transitions of the type a and b are performed, 
and after that we obtain non-vanishing 77 and hence non-trivial dependence on the new x 
variables. 

Also, adding or removing a polygon by transitions of the form of Fig. lg or h is 
straightforward. Because of the Hamiltonian and Lorentz constraints on the extra polygon 
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the wave function depends neither on the Lorentz orientation nor on the internal time 
variable of that polygon. 

Let us stress once again that it is the transitions that cause our system to be highly 
non-trivial. The classical system already has shown that infinite successions of such tran- 
sitions often occur (usually resulting in ever increasing values for pi , in which case the Li 
equally rapidly decrease. So the consequences of the constraints induced by the boundary 
conditions of Fig. 1 are severe. 

Ultimately, since we are performing cosmology rather than local quantum mechanics 
it is not so much the Schrodinger equation but rather a Wheeler-DeWitt type equation to 
which the entire wave function will obey. 

5. Discreteness 

The most striking consequences of the quantum structure of this model for some reason 
has never been observed or stressed by other authors. It is the discreteness of the relevant 
variables. First let us concentrate on the time variable. 

The Hamiltonian is the total deficit angle. For a closed S 2 universe this is An. Locally, 
the contributions to the Hamiltonian govern how parts of the universe evolve with respect 
to a time variable fixed at "distant polygons", or "distant observers". It seems that the 
very physical nature of our approach allows us to see this more clearly than otherwise. 
What we see is that the local Hamiltonians are also angular. Of course these angles are 
defined only modulo 2n , and so our Hamiltonians are also only defined modulo 2ti . Indeed, 
all expressions we have for the Hamiltonians in terms of the p\ give us only cos Hi (eq. 1.7), 
and to some extent also sin Hi (eq. 1.4). This means that what we really have is only direct 
expressions for e ±lH in terms of single- valued functions of pi. 

But this we consider as highly interesting. Apparently the evolution of the system 
is only well determined for integer time steps. Fractional time steps are ill-defined and 
skipped. Clearly time is quantized in our model. In fact, this quantization of time was 
seen earlier when the relation was established between angular momentum on the one 
hand and a time shift along a contour around a set of particles on the other. Since angular 
momentum is quantized, time shifts are quantized also [5]. 

Time quantization is also essential for a discussion of uniqueness and unitarity of the 
system. We want the evolution to be described unambiguously. One can only hope to 
obtain such an unambiguous law for time steps that are integer. Yet even there things 
are not quite this simple. Eq. (1.7) does give us the cosine of the contribution to the 
Hamiltonian, but not the sine. This means that, from that equation alone, we obtain the 
operator 
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a combination of a step to the future and a step to the past. What is needed is an 
unambiguous expression for e~ lH alone. 

Fortunately we also have (1.4). In the classical system this expression is sufficient to 
determine all angles uniquely at every transition. Given the rji there is in principle one 
overall sign ambiguity for the sines of the angles at each vertex. Of course what this means 
physically is that our quantum system evolves under laws that ere symmetric under time 
reversal. All that is needed is that for one diagram all angles must be given. We can then 
use eq. (1.4) to lift the sine ambiguity for all transitions to all other diagrams. 

To establish whether our system will be copletely unitary is still difficult however. 
If we take a complete cosmology we know that large series of transitions can relate one 
diagram to the same diagram at many different values of the parameters. So we cannot 
allow for all possible states ip(x) because most of them will violate the boundary conditions 
of Fig. 1, or, they will violate the Wheeler-DeWitt equation. It is at a local level that our 
equations seem to be sufficiently stringent to determine e~ lH completely. 

Since time is discrete one may now ask to what extent space is discrete also. It 
would be natural, perhaps, to suspect that space-time forms a complete lattice. But the 
situation is considerably more complicated. In principle the Xi can take any (integer or 
fractional) value. The thing to observe however is that the only p dependence comes from 
the hyperbolic sines and cosines of pi, or, we just see the operators e ±Pi occurring in our 
expressions both for the cosines and for the sines of Hi . But these are the operators 
for shifts over exactly a distance i , which is a unit step in the imaginary direction. The 
question now is how to exploit this fact. 

It is now of importance to use the constraint from (3.4) that pi are non-negative. This 
means that the wave functions will be analytic in the entire upper half plane, Im(xj) > . 

We have (we omit the indices j) 

*{z) = ±- r^^,iflm,>0. (5.2) 

^ V > 27TZ7_ 00 X-Z ^ ' 

Let us now define the amplitudes 

n = V(m) = -L l°° (n > 0) ; ^ = </>(0) . (5.3) 

2m j-oo x — in 

Then the operators e ±p act on these in a very simple way: 

e ±p (P n = (j> nTl . (5.4) 

We are now in a position to interpret the eqs. (1.4 — 8) as difference equations on the 
wave functions 4>{ nj }{{tF}) • There is a time variable at each polygon F , which we are 
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allowed to vary separately. In a given diagram each edge j has a (non-negative) variable 
rij . We lable the polygons and the edges in a cyclically symmetric way as in Fig. 2. There 
is now an algebra of operators. At each corner i there is an operator Tj = e~ lHi , defined 
such that 

T F = 11 Ti = e~ iHF (5.5) 

i€F 

is the time displacement operator at polygon F: 

2>0{n,-}(*l> • • • > tF, ■ ■ •) = 0{„,.}(*l, • • • ,t F + 1, • • •) • (5.6) 
And at each edge j there is an operator Uj = e Vi defined such that 

( -V »/;....(> 1- ■■■)= 0ni,...,n J —l...(*l> • • • ( 5 - 7 ) 

The entities Sj, Cj,cr,,-, 7^ in eqs (1.4 — 8) are now 



= -|i(Ti - I;" 1 ) , Cj = |(T, + T~ x ) , 



(5.8) 



Because of the observed analyticity in the upper half plane we have as a boundary 
condition 

0...,n-,... ^0 as rij — > +00 . (5.9) 

Now most of our equations will be of second degree (that is, involving at least two steps) 
in the rij direction, so that an other boundary condition may be needed. This is the 
value of the wave function if) at the origin. Here an edge vanishes, and hence the wave 
function must coincide with other wave functions for different diagrams, to be obtained 
via transitions as given in Fig. 1. This is the reason why in eq. (5.3) we considered only 
the set of wave functions that is connected to the origin by integer vertical steps. These 
are probably more essential than the ones we would have obtained had we started at an 
arbitrary other point in the complex plane. 

Finally we note that the equations at an edge to which a particle is connected, eqs 
(1.2) must be treated exactly as eqs (1.4 — 8). The mass m is an arbitrary free but fixed 
parameter here. Only one problem has not yet been addressed. This is the fact that, since 
the distances are now discrete, the distance between two particles can become exactly equal 
to zero. This was never a concern in the classical case, because such an event would occur 
with probability zero. Now it is a finite possibility. Presumably there will be room here 
to enter non-trivial non-gravitational interactions among the particles themselves, this we 
have not yet worked out. 
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6. Discussion and Conclusions 

Our procedure with tesselated Cauchy surfaces turned out to be strikingly suitable for 
a desciption of not only classical but also quantized particles gravitating in 2+1 dimensions. 
The lengths Lj of the edges of the polygons and the Lorentz boosts 2r]j across these edges 
turned out to be each other's canonically associated degrees of freedom, and the Poisson 
brackets, eq. (2.7), are quite suitable for setting up a quanization procedure. 

But replacing Poisson brackets with commutators must be done with care. Often, if a 
Hamiltonian is not quadratic in the momenta, a theory may turn out to become non-local, 
non-unitary or lacking a stable vacuum state. In this model we faced all these dangers. 
Now there seems to be a general rule that if a model is classically integrable it will have an 
integrable quantum version as well. It seems that this rule works to our advantage here. 

The effect of the quantization procedure is remarkable. Because the Hamiltonian is 
an angular variable the time coordinate comes out automatically being discrete. Only 
over integer time intervals the wave function evolves unambiguously. Because of Lorentz 
invariance one could have expected therefore that also the spacelike dimensions should be 
discrete. Instead, it is the imaginary parts of spacelike distances which will be quantized 
in integers. We found that the independent variables tp and nj in the wave function 
4>{n,j}({tF}) can all be restricted to integer values, because the wave equations in terms of 
these variables turned into difference equations. 

We have here an unusual analogue of the so-called Wick rotation in quantum field 
theory: instead of replacing time by an imaginary quantity we have kept time real but 
replaced the spacelike coordinates by imaginary numbers. Of course in both cases space- 
time obtains a (locally) Euclidean signature. 

Whether the variables Xj can be completely replaced by the rij remains to be seen. 
This replacement appears to become inefficient at large distances. A study of the possible 
classes of analytic wave functions reveils that there exist indeed functions that are zero on 
z = +in, n > and finite when z is real. These are necessarily divergent on the lower half 
of the complex plane. Therefore the functions <p are n °t completely representative for all 
states. One cannot build a complete basis out of them. 

A complete formulation of "quantum cosmology" in 2+1 dimensions has not yet been 
given. What we would like to see is an S'-matrix construction: given some asymptotic 
states \ip)in at time t — > — oo (if the universe started being infinitely large) or t = to (if the 
universe started with a big bang at t = t ), and out at time t — > +oo (for an expanding 
universe) or t = t\ (for a universe with a "big crunch" at t = ti), we would like to be able 
to compute the "scattering matrix" 

out(V#}in • (6.1) 

A problem here is to give a semiclassical description of the asymptotic states. This seems 
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to be all right if the universe became infinitely large there, but in the crunching case this is 
very problematic. In an earlier paper we expressed the suspicion that the crunching states 
become semiclassical also. This however was based upon the hope that the momenta as- 
sociated with the lengths Lj were something like the hyperbolic sines or cosines of the r] 
variables, which was not so crazy because the classical expressions for the momenta do 
contain the Lorentz 7 parameters. Now we know that this is different in the strong gravi- 
tational case, and the asymptotic states will keep their fundamental quantum mechanical 
nature. 

We actually found from the triangle equations (1.4 — 10) that 



is strictly bounded, even as 77 — > 00. One can conclude from that that during a classical 
crunch the quantity 



approaches a finite limit. Since this quantity counts the number of wave nodes one can 
deduce that the asymptotic crunching state cannot be described semiclassically. Therefore 
it will be very difficult to even define the matrix (6.1). We do not know how to characterize 
complete sets |V')in,out without overcounting. 

From the fact that (6.2) approaches a finite limit it also follws that the critical coeffi- 
cient k mentioned in Ref. [7] is equal to one. 

An alternative approach to a physical interpretation of the quantum theory could be 
to concentrate on the definition of an S'-matrix for scattering in an open universe, as was 
done in Ref. [5]. In that case the asymptotic states are always expanding, but the total 
energy must be constrained to be less than 2tc. Such a constraint would make it harder 
however to formulate the usual conditions of unitarity and causality for the S'-matrix. 

Discreteness of time has the consequence that energy is only defined modulo 2tt. We 
could call this unit the Planck energy. In other theories with discretized time this is a 
problem, because then there may not be a well-defined stable vacuum state [6]. Here we 
are not so much concerned with that. The total energy of the universe is only 47r, or two 
Planck units. So if we take a small section of this universe then the energy quantum is 
much too large to cause us concern. 

Can one add non-gravitational interactions to the model? What about a X(p 4 term 
in the Lagrangian? We must observe firstly that, although we seem to have a completely 
quantized model here, we have not yet seen creation or annihilation of particles. We know 
that creation and annihilation do occur in a non-gravitaional theory with a A</> 4 interaction. 
One should suspect this still to happen if one then adds gravitation. Secondly, so far 



3 




(6.2) 




(6.3) 
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we ignored the states where two particles coincide. Just because of the discreteness of 
the distances Lj this may be a serious omission that has yet to be addressed. Notice 
furthermore that the relevant equation linking energy and momentum, eq. (1.2c), only 
contains cos m , not sin m . So one may easily generate a sign difficulty for the particle 
mass m, comparable to difficulties that led to the necessity of second quantization in non- 
gravitational field theories. Since the transitions of Fig. lb, c and e leave no ambiguity for 
the deficit angle corresponding to the Hamiltonian we do not expect particle creation or 
annihilation to occur in pure gravity. But as soon as other interactions are incluided we 
probably will have to deal with a complete Fock space. 
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